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ABSTRACT 

We present analytic formulae for the integral number count distribution of 
cosmological bursting or steady sources valid over the entire range of fluxes, including 
density evolution and either standard candle or a power law luminosity function. These 
are used to derive analytic formulae for the mean redshift, the time dilations and 
the dispersion of these quantities for sources within a given flux range for Friedmann 
models with i7 = 1, A = without K-corrections, and we discuss the extension to cases 
with i7 < 1 and inclusion of K-corrections. Applications to the spatial distribution of 
cosmological gamma ray burst sources are discussed, both with and without an intrinsic 
energy stretching of the burst time profiles, and the implied ranges of redshift z are 
considered for a specific time dilation signal value. The simultaneous consideration 
of time dilation information and of fits of the number distribution versus peak fiux 
breaks the degeneracy inherent in the latter alone, allowing a unique determination 
of the density evolution index and the characteristic luminosity of the sources. For a 
reported time dilation signal of 2.25 and neglecting [including] energy stretching we 
find that the proper density should evolve more steeply with redshift than comoving 
constant, and the redshifts of the dimmest sources with stretching would be very large. 
However, the expected statistical dispersion in the redshifts is large, especially for 
power law luminosity functions, and remains compatible with that of distant quasars. 
For smaller time dilation values of 1.75 and 1.35 the redshifts are more compatible with 
conventional ideas about galaxy formation, and the evolution is closer to a comoving 
constant or a slower evolution. More generally, we have considered a wide range of 
possible measured time dilation ratios and discuss the values of the density evolution 
and the redshifts that would be expected for different values of the energy stretching. 
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Subject headings: cosmology — source counts — gamma-rays: bursts. 

1. Introduction 

The integral distribution N( > F) of the number N of sources with flux greater than F provides 
valuable information about the luminosity and spatial distribution of unidentified astronomical 
sources, especially in the absence of independent distance indicators. The main problem is that 
it convolves the information about the luminosity, density and distance distributions in a manner 
which is difficult to untangle. For gamma ray burst (GRB) sources, the angular distribution is 
highly isotropic and appears consistent with either a cosmological or an extended galactic halo 
interpretation (Fishman, et al. , f994). In either case, the departure of the integral distribution 
from a simple Euclidean N{> F) oc F"-^/^ law observed at low fiuxes may be an indication that 
the end of the spatial distribution has been reached (Meegan, et al. , 1992), and/or that at low 
fiuxes the effects associated with the luminosity function or the density evolution start to dominate 
the integral distribution (e.g. Wasserman, 1992, Wijers & Paczyhski , 1994). In the cosmological 
case, even an unbounded and unevolving standard candle distribution will slowly turn over at 
low F due to cosmological redshift effects z unity, e.g. Mao & Paczyhski , 1992, Dermer, 
1992. However, such effects at low fiuxes depend on the type of luminosity function and density 
evolution of the sources. 

Information about the redshift of the sources would be of significant interest for a cosmological 
distribution, as it could fix one of the crucial quantities which is otherwise an unknown parameter 
in statistical N{F) vs. F fits (e.g. Loredo and Wasserman, 1992, Fenimore, et al. , 1993, Band, 
1994, Horack, et al. , 1994, Emslie and Horack, 1994, Cohen and Piran, 1995, Meszaros & 
Meszaros , 1995, hereafter MM95). In the absence of identified counterparts and/or of identifiable 
lines, some information on the redshift could be obtained from the detection of a "cosmological 
signature", fainter (distant) sources being expected to have longer characteristic timescales due 
to cosmological time dilation (Paczyhski , 1992, Piran, 1993). Evidence for this effect has been 
reported from BATSE 2B data (Norris, et al. , 1994, Norris, et al. , 1995, Fenimore & Bloom, 
1995) but is not seen in some other analyses (e.g. Mitrofanov, et al. , 1994, Mitrofanov, et al. 
, 1995). One difficulty is that this effect is cleanest for standard candle sources with a standard 
duration; a broad luminosity function and/or an intrinsic spread in the durations could smear out 
the signature. Another possible difficulty with this signature is that its effects could be mimicked 
by intrinsic properties of the sources (e.g. Meszaros & Rees, 1993, Brainerd, 1994, Band, 1994, Yi, 
1994). An additional complication is that an intrinsic energy stretching of the time profiles may 
be present (Fenimore & Bloom, 1995, Fenimore, et al. , 1995) which would also weaken the 
cosmological signature. 

In this paper we present new exact analytic solutions for the integral distribution N{> F) 
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with a power law or standard candle luminosity function and a power density evolution law, in the 
context of cosmological distributions. These expressions are exact over the entire range of fluxes 
F, and thus are more general than the asymptotic forms discussed in MM95. We discuss mainly 
the i7 = 1, A = case without K-corrections (where i7 is the present ratio of the density to the 
critical one, A is the cosmological constant), with comments on extensions to the cases i7 < 1 
and inclusion of K-corrections. We also discuss analytic expressions for the theoretically expected 
redshift and time dilation, calculating the mean values and the dispersion under the effect of both 
a density evolution and a luminosity function. The cosmological signature is sensitive to both, 
the main effect of a luminosity function being a reduction of the signature over what is expected 
in a standard candle model. Some analytic expressions are valid both for bursting or steady 
sources, with an appropriate change of the evolution index. The inclusion of an intrinsic energy 
stretching of time profiles in the example of GRB (such that bursts are intrinsically narrower at 
larger energies) is discussed and incorporated into an analysis of the redshifts and cosmological 
signatures of GRB. We model the expected signatures both with and without energy stretching 
effects for the evolution and luminosity function parameters, and assuming a given value for the 
cosmological signature illustrate how one can determine the density evolution index and deduce 
the range of redshifts for sources in different peak fiux ranges. We also discuss the statistical errors 
associated with such redshift determinations, and discuss their compatibility with currently held 
views on the redshift of earliest galaxy formation. 



2. Integral Distribution 



In this section we derive analytic expressions for N(> F) of bursting or steady sources. 
Note that the corresponding differential distribution used for statistical fits is obtainable as 
N(F) =1 (dN(> F)/dF) \. The analytical relations collected in this section were already used 
(without equation details) in the fits to the 2B catalogue in Horvath, et al. , 1995 (hereafter 
HMM95). Note that our calculations in §2 and §3 are general, and hold for any cosmological 
sources with F in units of photons/(cm^s). Note also that in this paper we denote the fiux by F 
(unlike in MM95 and HMM95). This more general notation is motivated by two things: first, our 
calculations in §2 and §3 may be applied in some cases to steady sources; and second, even for 
gamma ray bursts the expressions are valid either for average or peak photons fiuxes, as long as F 
is in units of photons/(cm^s). We prefer F rather than the P sometimes used for peak fiux, which 
can be confused with power or period, and this also distinguishes it from the instrumental counts/s, 
which are denoted generally by C . Thus, in §2 and §3 F denotes generally photon number fiux, 
in the remaining sections photon number peak fiux of GRB, and C are the instrumental counts/s. 
It is also necessary to emphasize that our value of F is not to be confused with the energy fiux 
measured e.g., in units ergs/(cm^s). This is because the relation used here between the fiux and 
luminosity (see below equation (3)) holds only in the case when the luminosity has the dimension 
photons/s, and the fiux has dimension photons/(cm^s). Otherwise, for energy fiuxes our equation 
(3), and hence the following calculations, should be correspondingly changed (MM95). 
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We assume a density evolution as a power law of the scale factor, 

n{z) = no{l + zf , (1) 

where z is redshift, _D is a real number characterizing the density evolution, n{z) is the source 
proper density, Ug is the corresponding density at 2; = (-D = 3 corresponds to a constant 
comoving density), and the units of Ug are in Mpc~"^yr~^ for bursting sources or Mpc""^ for steady 
sources. The luminosity function $ is taken to be either of the standard candle type or a power 
law in the luminosity between some lower and upper limits in the photon luminosity L (in photon 
units s~"^), 

noS(C — Co) 1 (standard candle) ; 



^'^^^ 1 nC,^(CICm) , for £m < £ < Cm (power law) 

where n = Uoifl — !)/[! — for /3 7^ 1 and K = Cm I C,m- For i7 = 1, if there is no 

K-correction (MM95), the flux F is related to C and the comoving radial coordinate x (Weinberg, 
1972, MM95) by 

F = Fh{1-x)\-' , FH = C/{47rRl) 
1 + z = (l-x)-' = v~' , (3) 

where x = ^ — 1] varies between and 1 for z ranging from to 00 (77 is the conformal time), and 
where Ro = (2c/ Ho) = 6000/i~^Mpc is the Hubble radius (c is the velocity of light; Ho = lOO/i 
km/(s Mpc) is the Hubble parameter). 

For a given density evolution index D the number of observed sources with flux greater than 
F (the integral distribution) is given by (MM95) 

Nn{> F) = AirRl / ^C)dC / (1 - x)'^''' dx , (4) 

■JCm -JO 

where for bursting (steady) sources B = 1 (B = 0), 

Xr = il + [47rRlF/C]'/y = (1 + a^/^ , (5) 

and we defined the normalized peak fiux a = F/Ffj as the ratio of the observed fiux to the 
Euclidean fiux from the same source at the Hubble radius. Concerning the parameter B, this 
is introduced in order to treat both steady (B = 0) and bursting sources (B = 1). The latter 
value arises because in this case one must include an extra (1 + z)~^ factor into equation (4) to 
account for time dilation of the burst rate per unit time, e.g. MM95, Mao & Paczyhski , 1992. In 
what follows, we will consider bursting sources, i.e. B = 1. Nevertheless, the formulas presented 
in the entire §2 and §3.1 apply also for steady sources; for them one should use the appropriate 
dimension for steady sources and in §2 and §3.1 use the results based on equation (4) with B = 1, 
but substituting _D by _D + 1. This means that what holds for a given D in bursting sources also 
holds for (D — 1) in steady sources. In both cases, D = 3 means comoving constant source density. 
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Analytic solutions of equation (4) were given by MM95 for F in the two (three) asymptotic 
regimes of the standard candle (power law) luminosity functions (2) and density evolution (1) in 
the special case of _D < 4 with D integer or half-integer. It is useful, however, to have analytic 
solutions which are valid for arbitrary F and D. In the next two subsections we derive the full 
analytic solutions over the entire range of F and D. 

2.1. Standard Candle Luminosity Distribution 

For the standard candle (SC) case in equation (4) we have ^(C)dC = Kq, so the 
double-integral reduces to an integral over the single variable x- In MM95, the solution of the SC 
case was found by expanding the binomial in the integrand for integer or semi-integer values of 
the density evolution index D < 4. Here we do the direct integration of equation (4) in the SC 
case using r] as the variable of integration. The result is 

Nd(> F) = (47r/3)n„i2>;3/2 j^^^^^ ^ 

-3/2 . ^ /l-^?"''' l-r]l'-^''^ 

^^^^"^ = ^y^^ 5^r^+ 11-2D ) ' 

Yim^Noiy F) = (4:TT/3)noRlAD , D < 4.5 

lim Nd{> F) = (47r/3)n„i2>;3/2 ^ 

Ad = Q/[(9-2D)(10-2D)(ll-2D)] , (6) 
where Id is a dimensionless function of (Tq or r]i which are defined as 

(To = F/Fho = F/[£o/{47rRl)] ; r^, = 1 - xi = [l + a;'^^] (7) 

The expression of N^iy F) in (6) is of broader applicability than that in MM95, being valid 
for any real D (except D = 4.5, 5, 5.5). (In MM95 only the two limiting cases were obtained for 
arbitrary D < 4, and the identical results in different forms for integer and semi-integer D < 4.) 



For the three remaining values of D the direct integration of equation (4) gives 

(t;'/'/4.5 = -(9/2) + 6771 -3 In - (3/2)772 , (8) 
(T;3/2jg = 377-1 + 6 In 771- 3771 , (9) 
(t;'/'/5.5 = (9/2) - 3 In 771 - 677fi + (3/2)77^2 . (10) 

Note that for D > 4.5 one has limj?^o A^d(> F) = 00. Note also that for _D = 6 we have 

Ne{> F) = {47r/3)noRl[xi/{l - Xi)f = {47r /3)noRla;^^^ , (11) 



i.e., /g = 1. This is exactly the same expression as in the Euclidean limit for arbitrary D (see 
equation (6)). For D < 6 (D > 6) the expressions N^iy F) grow less (more) steeply towards 
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small F than the corresponding Euclidean curve (11). In other words, < 1 for D < 6, and 
Id > 1 for D > 6. For example, for D = 7 we have 

Nr(> F) = (47r/3)n„i2>;3/2 [l + (3/2)a;'/' + (3/5)a;'] . (12) 

Note that the mimicking of an Euclidean behavior for an evolution law _D = 6 is characteristic 
of an i7 = 1, A = bursting model. For instance, with the same model but steady sources, 
the Euclidean behavior occurs with D = 5, rather than 6. More generally, for any arbitrary 
cosmological model one can always find an evolution law that will just cancel out the cosmological 
effects. 

In Figure 1 we show two sets of theoretical curves for two values of Cg- For comparison, we 
also show the 2B observed integral numbers, with a sliding vertical axis. The top set of curves 
has an Cg chosen to give an approximate eye fit of the 2B data to the D = 4 curve, while the 
lower set has an Cg giving an approximate fit of the data to the D = 2 curve. While this is not an 
accurate fit, it illustrates the fact that fits to observed data can be found for various D (including 
_D = 3, not shown) by varying Cg (and ng, but the vertical axis has been left arbitrary). More 
detailed joint fits to the 2B and PVO data sets give somewhat different D,Cg best fit values 
(see IIMM95, and Table 1). Fits to 3B data are in progress. Qualitatively, this D degeneracy can 
be understood by considering, e.g. what happens as one increases Cg, in which case one needs to 
see deeper to maintain a given fiux F. Seeing deeper gives more cosmological bending to N(> F) 
for nonevolving bursts; to match the data one would then require an N(> F) that bends less 
severely by incorporating a density evolution that increases the relative frequency of bursts at 
large distances. 



2.2. Power Law Luminosity Distribution 



The equation (4) for a power law (PL) luminosity function given by the second line of 
equation (2) reduces to the form (see MM95, equations (14)-(17)) 

NuO F) = In = fnR^ga-Jf^ I, , (13) 

where 

„6A'i/2 8-2D 

Ij, = 2b-'+'^ / dy y'-'^{l + y)-'[l + ^ a^y^l + y)-'] , (14) 
•^^ k=i 

and Id is again a dimensionless function of a^, K , [i and _D, valid for integer and semi-integer 

values of _D < 4 (for other values of D see below). We have defined here 

y = b(CIC^fl\ b = a-'/^ = [£^/{4Tr RlF)]'/^ a^ = FlFHm 

Fum = Cm/i^TrRl) , FhM = CM/i^TrRl), K = CM/Cm > 1, 

ak = (-l f[3/(k + 3)][(8-2Dy./[kl(8-2D-ky.]l A; = 0, 1, (8 - 2i?) . (15) 
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No is a function of F ((Jm), and depends on the parameters ra, ^M, P and _D. In MM95 

the integral (14) was estimated in three asymptotic regimes of F for arbitrary (3. However, exact 
solutions valid over the entire range of F are of interest and in principle can be derived for any 
rational value of (3. An exact solution is particularly simple for some definite rational values of 
(3. Guided by the value (3 ~ 1.88, which fits the slope of the low F differential distribution (e.g. 
Meegan, et al. , 1992), we adopt here [3 = 15/8. A nontrivial choice of [3 must be near such a value, 
where the shape of the distribution at the faint end is dominated by the shape of the luminosity 
function, rather than cosmological effects. On the other hand, significantly steeper (e.g. (3 2) 
power laws are dominated by the faint end, while significantly fiatter power laws (e.g. (3^1) are 
dominated by the bright end, and thus mimic a standard candle behavior, the faint part of the 
distribution being dominated by cosmological effects. With this value of /3, the expression for lo 
in the integer and semi-integer case D < A becomes 

/^ = 85-^/M dtt^F-^J2^k{tVF) =8b-'^^AD^M. (16) 

where 

t = y^/^ = b'/\C/C^f/\ F=l + t\ (17) 
= b'/^ = a-'/\ tM = b'/'K'/' = iFHM/Ff/'; 
Ad = 6/[{9-2D){10-2D){ll-2D)] , M = 8 - 2D , 

and (see Appendix) 

ftM ^ ftM ^ 

^M = A^' dtt*F-'J2 ""kdVFf = dtt*F-'J2i(k + l)(k + 2)/2]F-K (18) 

im J n im J n 



The $M have an exact solution over the entire range of F, and can be obtained recursively starting 
from $0 and $i as shown in Appendix. The result for the indefinite integrals is 

Mt) = (3/8)(l/16v^)0(O + (l/32)(i/r)-(l/8)(i/r2) , 

Mt) = (l/4)(iVr') + (33/12)V'o(0 , (19) 

-2 



where the definite integrals are $m = i^Mi^M) — i^Mi^ra)-) ^iid 

^2 _|_ ,yf^j^ _|_ 

Q(t) = In ^ + 2 arctan(V2i + 1) + 2 arctan(y2i - 1) . (20) 

To get the asymptotic regimes, note that the function AD'tpMi't) = I dt t'^F~'^ Ylh=o O'kit'^lT)^ 
for < i ^ 1 may also be written as a power series of t. Its concrete form may be obtained by 
first using the formula for the sum of geometrical series, as well as F~^ = X]s^o(~l)*^^*5 and 
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then doing the integration. This gives Aui'Mi't) = (^^/5) + ••• (the integration constant is zero, 
because t/'M(0) = 0), and for any D the first term is the same. Then, if tM ^ 1 is so small that 
for t = tM one may restrict oneself to this first term, one immediately obtains the Euclidean limit 
with Au^M — Hence, one reproduces the asymptotic t ^ Q {F ^ oo) behavior obtained 

in MM95 (see equation (18) in MM95). Similarly, for t 1 the function Ao'i'Mii) may also be 
written as a series, if one uses = t~'^{l + = t~'^ Z]^o(~l)*^~^*- Then, after integration, 
one has A£)ipM{t) = Auiconst. — (7t'^)~^ + ...). (Here the integration constant is not zero, and 
is identical to lim^^oo V'm(05 which is easily calculable using lim^^oo a'rctan(-v/2^ ± 1) = ir/2. 
Nevertheless, this integration constant is cancelled, once the definite integral $m is calculated.) 
H ^ 1 is so big that for t = t.^ one may restrict oneself to this first term, one reproduces the 
second asymptotic regime of MM95 for t oo, F Q (equation (24) of MM95), and one has 
Nd(> 0) = (Att /SjnoRoAu, as it must be. The third asymptotic behavior predicted in MM95 
occurs for values of K large enough that in some range of F one may simultaneously approximate 
t = tm with the Euclidean limit, and t = tM with this second limit (see equation (20) of MM95). 
Thus, for large K the formulae in this section also give the three asymptotic regimes described in 
MM95; but for K ~ 1 the situation is similar to the SC case, as expected. 

For the case _D > 4, when D is integer or semi-integer, in equation (13) Id are also directly 
calculable. First, one must solve the x integrals in equation (4) directly, as in the SC case, and 
then one integrates over C The results are 

I4.5 = h-^'l^ [(24/7)^-3 - (24/7)i-^lnr - MF-^ + (9/[28v^])0(Ol , 

h = 5-'/^ [-(48/7)^-3 + (48/7)i-^lnr + (6/[7V2])0(Ol*'' , 
h.5 = 5-'/Ml2i + (24/7)^-3 - (24/7)i-^lnr - (24/[7v^])0(d*'' , 

le = (8/5)(A'^/«-l), (21) 
/e.5 = (8/5)(A'^/«-l) + (3/4)(A'9/8_ 1)^-1/2^ 

Ir = (8/5)(A'^/« - 1) + (4/3)(A'9/8 _ 1)^-1/2 + (24/Q5)(K''/^ - l)a-' . 

We see that Iq does not depend on dm, and hence Nq(> F) oc oc F~^^^. Thus for D = Q 
we obtain the Euclidean case both for SC and PL luminosity function. In fact, one can easily 
show that for D = 6 one obtains the Euclidean case for arbitrary luminosity functions ^(C). From 
equations (4) and (5) it follows for D = 6 

Ne(> F) = (47r/3)(47rF)-3/2 / ^(jC)jC^/^dL = const x F'^^^ (22) 

Jo 

where const does not depend on F. Obviously, for any ^(C) and D > 6 the expected integral 
distribution N£)i> F) is even larger than the Euclidean value. As mentioned at the end of §2.1, the 
Euclidean behavior of _D = 6 is for i7 = 1, A = and bursting sources. For arbitrary cosmological 
models, it is always possible to find an evolution law that cancels out the cosmological effects. 
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Note that in Section 2.2 we restricted ourselves to the integer and semi-integer values of D. 
Nevertheless, the analytical formulae of this Section may be generalized for any real D. The 
relations (4-5) hold for any real D. Then the key problem is to write down the definite integral 

/ X^il-xf-^''dx= / QD{x)dx. (23) 
Jo Jo 

The primitive function of Qd{x) foi" non-integer (8 — 2D) is obtainable in at least two different 
ways. First, one writes for M = 8 — 2D the Taylor series 

oo 

(1 - X)"^ = 1 + - - k + (24) 

k=i 

where k is integer. This is not new for M integer (see MM95), where the series is finite because 
for k > M -\- 1 one obtains identically vanishing terms. On the other hand, for non-integer M we 
have an infinite series. In any case, one should use the relation (5), and in this series one obtains 
integer exponents of £. Then the integration over L is analytical (although it can be cumbersome) 
and in principle can be done for any rational [i. Second, one substitutes (1 — x) = 77, and obtains a 
simple primitive function of Qo from integration over 77; this primitive function is a sum of three 
components, containing either In 771 or - in general non-integer - powers of 771. Therefore here, for 
non-integer M, the technical problems can arise from the integration over L. In other words, the 
simpler integration over x may lead to a more complicated integration over L. In this paper we 
do not discuss these technical questions any further, and restrict ourselves to the cases of integer 
M = ^-2D. 

In Figure 2 we show two sets of theoretical PL curves, both for K = Cm / = 10^. The 
particular choices of CM^^m used here are not unique, but have been selected so as to give an 
approximate eye fit of the 2B data set to the D = 4 curve (top group of curves), and to the 
D = 2 curve (lower group of curves). The visual match of these theoretical PL curves to the 2B 
data is as good as that for the SC case in Figure 1, even though here we have a luminosity ratio 
of 100. (In fact, detailed fits to joint 2B and PVO data sets with PL models of K = 10^ are 
indistinguishable to within la from corresponding SC fits, as shown in IIMM95 (see also Hakkila, 
et al. , 1994); such detailed PL fits to both 2B and PVO lead to more accurate choices of Lm^ D, K 
values, which have been used in Table 3). The approximate eye fits to 2B shown in Figure 2 
merely illustrate the fact that in the PL case, different choices of CM^C^m (even for constant K) 
lead to a different D which matches the same observed data set. 

In Figure 3 we compare the generic shape of the SC curves for various D (top set of curves) 
to that of the generic PL curves (lower set of curves). For the latter we have chosen an artificially 
large spread of luminosities K = 10^ with maximum luminosity Cm equal to that of the SC 
luminosity Cg- This large spread K permits one to see more easily the three limiting asymptotic 
regimes of the PL curves, going from the Euclidean F"-^/^ to F^~^ (where (3 = 15/8 is used) 
to _F ^ constant for _D < 4 as _F is decreased. The intermediate regime present in the PL case 
(e.g. Wasserman, 1992, MM95) is not present in the SC case. We note that the detailed x^ fits of 
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HMM95 indicate that there are statistically acceptable fits to the 2B data set alone with K = 10^, 
but only to K ^ 10^ if 2B plus PVO is used. 



3. Redshifts and Time Dilation 

3.1. Mean Redshifts and Dispersions 

The redshift of a source of luminosity C producing a fiux F is, from (3), (5), 

1 + ^ = [1 + (F/F^)-i/2]2 . Fh = jC/iAirRl] . (25) 



For a given ^(C) one has 

r^M rxi 

'Cm Jo 



, f^M fXl 

Nd{> i^i) = 4TrRl / ^£)d£ / QD{x)dx , (26) 
JCm Jo 



where Qd{x) — ~ x)^~^^X^ ^"^^ foi" the fiux Fi we have = [1 + (-Fi/-F//)^''^]~^. The number 
of observed sources, that give fiuxes between Fi and F2 is (0 < < < 00, including Fi = Q 
and F2 = 00) 



Nd{Fi,F2) = AttRI / $(£)rf£ / QD{x)dx = Nd{> i^i) - Nd{> F^) , (27) 

J Cm JyI 



'■Cm rX2 

'Cm Jxi 

where = [1 + (-^2/-?// )^''^]~^- In this equation there is a double-integral of the function 
^(C)QDix) depending on variables x ^- If takes the integral of the function 
(1 — y^)~'^^{C)QB{x) the same area of the variables £ and x, and divides the result by 
N£)(Fi, F2), obviously one obtains the average value of the quantity (1 — x^)~^- But since 
(1 -|- 2;) = (1 — x)~^, this gives the average value of the scale factor (1 + z) for all the sources with 
fiuxes between Fi and F2, namely 

, ltHndi:j-^{l-x)-'QD{x)dx Nn^,iF,,F2) 

^ ^ ltHC)dcI-:Quix)dx ~ ^M.F2) ■ ^''^ 

This expression follows from the relation (1 — x)~^QDix) = QD+iix) is valid for arbitrary 
luminosity functions. 

Also, since Qd+2(x) = (1 ~ x)~'*Q-d(x) (1 + — ~ x)""*? it follows similarly that 

n|,(PP)). /l-<i>(/:)^/:/-(i-x)-^Qz.(x)^x wf,,^,) 

^ ^ ' - !tHi:)di:i-:Qn{x)dx " mF,,F2) ' ^''^ 

Then 

A(l + z{Fr,F2)) = [(1 + z{Fr,F2)y - (1 + ^(^1,^2))'! (30) 
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characterizes the dispersion of (1 + 2;) around the mean value. Again, it is worth mentioning that 
these specific expressions are valid for i7 = 1, A = models without K-corrections. However, 
similar expressions can be derived for arbitrary models. 

Clearly, for a given Fi and F2 the values of (1 + z(Fi,F2)) and (1 + z(Fi, F2)y are different 
for different D. In the SC case they depend also on £0, and in the PL case on and K , but 
they do not depend on Kq, because the latter parameter cancels out in the definitions. 

As a consistency test of these relations consider a special case where F2 = (Fi + dF), where 
dF is infinitesimally small, and - in addition - we have an SC luminosity function. In this special 
case it follows that 

Nd(Fi, (Fi + dF)) = I (dND(> F)/dF) \f=f^ dF, (31) 
l + z(FiAFi + dFj) = (1-xi)-'; A(l + z(FiAFi + dF))) = 0, 

where we used the relation (dN£)(> F)/dF) = (dN£)(> F) j dx){dxl dF). We see that in this 
special case, the mean redshift equals to the redshift defined by xi = -\- (Fi/ Fh)^^^)~^ , and 
there is no dispersion, as expected. 

3.2. Time Dilations and Energy Stretching 

While the previous sections were valid both for steady or bursting sources, from here on we 
specialize to the latter case. Another observable quantity of potential cosmological significance 
for bursting sources is the intrinsic duration of source events, e.g. the time interval between the 
beginning and end of a burst, or a subpulse of a burst, or a characteristic variation timescale of 
the source fiux (or the inverse frequency of a spectral line, if such exists). This duration may be 
quite different for different energy bands. 

In the simplest case the intrinsic duration is "grey", or independent of the energy band, e.g. if 
there is a "standard" GRB burst duration Atg in the rest-frame of the source. Then the apparent 
duration of a source placed at redshift z is At = Ato(l + z) (Paczyhski , 1992, Piran, 1993), and 
two bursts placed at different redshifts zi and Z2 will show different apparent durations related by 
(Fenimore & Bloom, 1995) 

ri2 = iAh/At2) = (1 + ^i)/(l + Z2) . (32) 

A less simple case is that where the intrinsic burst durations are not grey, i.e. there is an 
intrinsic energy stretching, such that the intrinsic duration of the event in its own rest frame 
depends on the photon energy or waveband. This is reported to be the case with GRB, the 
same burst exhibiting increasingly shorter durations at higher photon energies (Fenimore, et al. 
, 1995, Fenimore & Bloom, 1995). As argued by these authors, the BATSE 2B data indicates 
an energy dependence of the type Ato(E) ~ Ato(Eo)(E / Eo)~^ , where Eq is an arbitrary energy. 
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Ato(Eo) has the dimension of time, and k ~ 2/5. Assuming this to be the case and - without loss 
of generality - taking k to be an arbitrary real number, one can show that if bursts had some 
standard duration at the same energy in their own rest frame, then the result of observing from 
different redshifts would lead to an apparent duration-redshift relation given by 

{Ah{Eo)/At2{Eo)) = [(1 + zr)Ato{Ei)]/[{l + Z2)Ato{E2)] = [(1 + ^i)/(l + ^2)]'"' , (33) 

where Ei = (1 + zi)Eo, E2 = Eo(l + 22). Because the right-hand-side does not depend on Eg, 
the ratio of apparent durations is [(1 -|- zi)/(l + Z2)Y~^ for arbitrary Eq^ and this energy need not 
be specified later. Thus in general we can write 

ri2 = iAh/At2) = [(1 + ^i)/(l + Z2)Y = , (34) 

where T12 is the time ratio, ri2 is the redshift factor ratio, j = 1 if there is no energy stretching, or 
J = 1 — A; if there is energy stretching. 

If we take events at two different fiux levels and Ff, (Ff, > Fd), assumed to come from an 
SC luminosity function, these correspond to redshifts z^, Zf, (z^ > Zf,) which come from the relation 

Fb,d = Co/[i7rRli^l + Zb,d - if] , (35) 

(valid for i7 = 1, A = 0, no K-correction; see §5 for a relaxation of some of these assumptions). 
Then, if we measure the ratio of the corresponding event durations (assumed to be some standard 
duration in their respective rest frame), T^b = (Atd/Atf,) > 1, we have 

(1 + Zd)/il + Zb) = tII\ Fb/Fd = [(TTT^ - l)/(7rT^ - 1)]2, (36) 

and hence 

VTTTb = - i)/(y^ - r]/^^^-)) , 

VTTVd = T'J^''^VTTTb . (37) 

Thus in this special case, both relevant redshifts are obtainable directly analytically. We see that 
one must have F^j Ed > t^/'' > 1, and as Eh/Ed runs from t^^'' to 00, zj, decreases from 00 to zero, 
i.e larger Eh/Ed leads to smaller zj,. On the other hand, for a fix Eh/Ed a smaller t^^'' gives smaller 
redshifts. 

To improve the statistics, one would want to bin the bursting source fiux into some 
appropriately defined ranges of "dim" and "bright" categories (taking only sources whose fiuxes 
are inside the ranges Ef, ± AFf, and Ed ± AEd centering around some mean values of dim and bright 
fiuxes Ed and Ff,). One can then use the formulas (28) derived above to connect the corresponding 
ratios of the mean durations (assumed to be standard, with or without intrinsic energy stretching) 
to the ratio of the mean redshifts. 



(1 + Zd) 



= (l + z((Ed-AFd)AEd + AEd))), 
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(1 + ^6) = (l + z((Fb-AFb)AFb + AFb))) , 

A(l + zd) = A(l + z((Fd - AFd) A Fd + AFd))) , 

A(l + zb) = A(l + z((Fb - AFb) A Fb + AFb))) , 

(l + zd)/(l + zb) = = (Kul^fl^ = tH' . (38) 

One can also calculate the dispersion in this last ratio of mean values by combining the relative 
dispersions (30) for both groups alone via 

A((l + ^d)/(l + ^6)) ^ r 2 } x2 

/ =-^ = J A(l + ^,)/(l + ^,) + A(l + ^,)/(l + ^,) . (39) 
(1 + Zd)l{l + Zb) rdb y ^ ^ ^ ' 

This can be done for different parameters and separately for SC and PL luminosity functions, the 
mean redshifts and dispersions being connected to the integral number expressions Nd for each 
type of luminosity function via equations (28,30). 



4. Application of Time Dilation Analysis 



Attempts to detect a "cosmological signature" (Norris, eA al. , 1994, Norris, eA al. , 1995) in 
GRB have compared the mean duration At of events in the BATSE data base with peak fluxes 
corresponding to count rates within specific bands C,id i AC'dd, Cd ± AC'd and C'b ± AC'b, where 
the subindices are "dd" for dimmest, "d" for dim and "b" for bright bursts. The values chosen by 
Norris, et al. , 1994 are (C'dd - AC'dd) = 1-4 Kct/s, (C'dd + AC'dd) = 2.4 Kct/s, (C'd - AC'd) = 2.4 
Kct/s, {Cd + ACd) = 4.5 Kct/s, {Cb - ACb) = 18 Kct/s, {Cb + ACb) = 250 Kct/s, where Kct/s 
= kilocounts/ second. These count values can be roughly translated (based on averages of data 
from Nemiroff, 1994, private communication) into the following peak fiuxes: 1.4 Kct/s 0.45 
photons/(cm^s); 2.4 Kct/s 0.70 photons/(cm^s); 4.5 Kct/s 1.15 photons/(cm^s); 18 Kct/s 
5.0 photons/(cm^s); 250 Kct/s 45.0 photons/(cm^s). In what follows we lump together the 
"dim" and "dimmest" category into a single "dim" one labeled by the index d' (standing for d+dd 
in the above categories). Hence we have 

Fd' - AFd' ~ 0.45 , Fd' + AFd' ~ 1.15 , Fb - AFb ~ 5.0 , Fb + AFb ~ 45.00 , 

Fd' ~ 0.80 , AFd' ~ 0.35 , Fb ~ 25.0 , AFb ~ 20.0 (40) 

in units of photons/(cm^s), implying a ratio Fbj F^t ~ 31. Note that these are approximate 2B 
bands, which may differ depending on the analysis and the data cuts, but they serve to illustrate 
the effects discussed. A more detailed discussion of these bands is presented in Horack, et al. , 
1995b, who use a ratio Fb/F^i ~ 21. 

In Norris, et al. , 1994, Norris, et al. , 1995 an observed ratio of durations Td'b — 2.25 was 
reported. More recent values are 1.75 (e.g. Norris, 1995b, Horack, et al. , 1995b) and 1.35 
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(Fenimore, 1995b). A first quick and simple estimate of the implications of such a time dilation 
can be obtained using the analytic formulae based on the fiux bin centroid values (37). This gives 
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(41) 

We see that the presence of an energy-stretching has an essential impact on the values of relevant 
redshifts; taken at face value a time dilation signal of 2.25 could imply that the dim burst are 
at redshifts larger than the most distant known objects (see also Fenimore, et al. , 1995). On 
the other hand, the lower time dilation values of 1.75, and especially 1.35, give more reasonable 
dim+dimmest redshifts, even with energy stretching. However, such a simple estimate based on 
equations (37) neglects binning effects as well as luminosity function and evolution effects. 

More reliable conclusions may be obtained by calculating (1 + z^t), (1 + 2:5), their ratio, and 
the corresponding dispersions via the formulas (28,30,38,39), which specifically include luminosity 
function and density evolution effects via the No entering the definitions of the these bin-averaged 
quantities. For a putative measured time dilation t^i}, one can then search for the cases when the 
ratio Td'h = (1 + Z(ii)l{l -\- Zf,) is equal (for j = 1, no stretching) to T^'h (2.25, 1.75, 1.35 in the 
Norris, et al. , 1995, Norris, 1995b, Fenimore, 1995b cases), or else is equal (for j = 3/5, with 

5 /3 

stretching) to tJ^^ (or 3.86, 2.54, 1.65 in the same three cases just mentioned). To do this we need 
to use the appropriate pairs of corresponding values Cg and D (for SC) or triplets of corresponding 
values Lm, K and D (for PL) which give good logiV — log_F fits (HMM95). We can restrict 
ourselves to values of _D < 5, since _D > 6 is excluded by the downward curvature of the counts 
at low F, and D = 5.5 gives also wrong fits (HMM95). The values calculated from our analytic 
expressions (38,39,28,30) using the SC values of No from §2.1 are shown in Table 1. 

These redshift estimates are systematically larger than the values from (41). Nevertheless, 
the dispersions of the average redshifts and their ratios are large, and therefore as a rough first 
approximation the values from (41) are compatible with those of Table 1. For instance, even in 
the case of t^i}, = 2.25, subtracting a la dispersion from the mean value of r^t}, for D = 5 gives an 
r^/f, comparable to the mean r^^f, value for D = 4.5, corresponding to (1 -|- Zd^) — 1 ~ 5.5. This is 
near the redshift of the most distant quasars known. 

The uncertainties in the mean redshifts are increased also by the fact that the estimates 
depend also on Co- The values used for each D are from fits to the BATSE 2B data by HMM95, 
and a ±1(7 variation around the best fit can change Cg by factors ~ 2. In Table 2 we illustrate the 
impact of such ^ la changes in Co for a fixed D = 4.5. We see that the impact on the redshifts of 
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varying between different allowed Co for the same D is non-negligible. Nevertheless, while keeping 
such uncertainties in mind, it seems clear that if the time dilation of 2.25 (Norris, eA al. , 1995) 
were correct, under the SC assumption Table 1 indicates 'zl ~ (0.2 — 0.6) and ~ (1.5 — 2.5) for 
j = 1 (no energy stretching), and ^ ~ (0.6 — 2.0) and 1^ ~ (4 — 10) for j = 3/5 (with energy 
stretching). For this time dilation, extremely large redshifts seem to be indicated for the dim 
objects (see also Fenimore, et al. , 1995, Horack, et al. , 1995b), but because of the dispersion 
around the mean values, even for Td^f, = 2.25 one cannot conclude necessarily that z^^ > 5 is 
required. For the lower values T^if, = 1.75, 1.35 the dim redshifts are of course more moderate. 
This is discussed further in §6. 

For the power law (PL) distribution luminosity function, we can use the same formalism 
to calculate the mean redshifts and their dispersions. We choose a relatively broad luminosity 
dispersion K = {LmI J^m) = 100, which while still giving a good fit to the 2B-|-PV0 data 
(HMM95), accentuates the differences between the SC and PL luminosity functions. While 
formally one uses the same equations (38, 39, 28,30), the values of No entering into them are now 
obtained from an integration over luminosities placed at different redshifts (§2.2). The results are 
summarized in Table 3. 

These values appear slightly larger than in the SC case, but no far reaching conclusions may 
be drawn from this difference. In fact, the remarkable difference between Table 3 and Table 1 
is provided by the much larger redshift dispersions for the power law luminosity distribution, as 
opposed to the standard candle case. Because of this large dispersion, the PL and SC results 
are clearly compatible with each other. Also, this very large dispersion implies that even for 
Tf^/ft = 2.25 the dimmest redshifts encompass within their la error bars the values z^i ~ 1, and the 
same true to a larger extent for the lower values of Td'b- The larger dispersion of mean redshifts in 
Table 3 is quite reasonable, because for SC at a given F one has a given z, but for PL at a given 
F there is range of z corresponding to a range of L. The implications of such large dispersions for 
cosmological models for a given measured time dilation are further discussed in §6. 

The dependence on in the PL case is illustrated in Table 4. We see that, similarly to the 
SC case, the dependence on is non-negligible, and while not changing the means by much it 
can change the ratios and further increase the dispersion. The dependence on K is illustrated 
in Table 5. Here D and are fixed, and K varies. In this case the r^'b practically remains 
unchanged, but the mean redshifts and dispersions vary significantly. 

5. K-Correction and i7 < 1 Effects 

In the previous sections the calculations of the redshifts of GRBs or other sources was done 
without any K-corrections and for i7 = 1. In this section the impact of relaxing these assumptions 
is estimated (c.f. also Fenimore & Bloom, 1995, Horack, et al. , 1995b). For simplicity, we 
illustrate these effects for the special case when equations (36, 37) hold, since it was found that 
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qualitatively they reproduce the behavior found in the more detailed analysis. We further simplify 
the situation by assuming that only one or the other effect operates alone, i.e., we either have 
J7 = 1 and a K-correction, or i7 < 1 and no K-correction. 

If £j, ~ z/"~^, where Lj^dv is the peak-flux luminosity in the frequency interval z/, (z/ + dv) 
and a is a real number, then there is no K-correction for a = 1 (e.g. MM95). If a 7^ 1 then the 
peak-flux F from a source at redshift z will be (1 -|- z)"~^ times greater than the expected peak 
flux when no K-correction occurs. (For instance, for a = and £j, ~ one will have a (1 -|- 2;) 
times smaller peak flux than for a = 1.) In other words, instead of equation (3) here one has: 

F = £/[47ri22(l + ^)i-«(yrT^- 1)2]. (42) 

Assuming two sources with the same intrinsic peak flux spectra characterized by a at redshifts Zf, 
and z^i , instead of equations (36, 37) one obtains 

(1 + ^d')/(i + ^b) = r]f , ^6(1 + z,f-"l[FAi + zd.f-"] = [(TTT^ - i)l{VTTT, - (43) 

and hence 

Vrn^ = ( V^^rir'^^V^.' - l)/( V^^rir'^^V^.' - rj^h VT+^' = rl!i'^^VT+T,. (44) 

The impact of a 7^ 1 can be seen from the fact that these equations are identical to equations (36, 
37), if one formally substitutes Fh/Fd' by Fj^t^^ /Fd'. Hence we see that for a < 1 (a > 1) 
one obtains systematically smaller (bigger) ratios for the peak fluxes, and hence bigger (smaller) 
redshifts than in the absence of a K-correction. To illustrate this effect consider again Ff, = 25 and 
F(ii = 0.8 (in units of photons/(cm2s) ). For a = 1 we had (equation (41)) Zf, = 0.6 and z^i = 5.2 
for j = 0.6. For the same case, but with a = (i.e. with ~ '^~'^)^ we obtain Zf, = 1.1 and 
Z(i' = 7.2. This means that one may expect that K-correction effects, when relevant, will increase 
the redshifts that one derives. Generally a = 1 is a good fit for the peak fiux spectra of GRB in 
the range where most photons are collected from (e.g. Band, 1994), and therefore K-corrections 
may not play a large role for most sources. Nevertheless, one cannot exclude the situation where 
a systematically smaller a may be needed (e.g. Schaefer, et al. , 1994) where K-corrections are 
relevant . 

For z ^ few a value of i7 7^ 1 may also affect the estimates of the relevant redshifts. We 
restrict ourselves to i7 < 1. For a = 1 we will have for a source with peak luminosity C at redshift 
z a peak fiux 

F = jC/[4Tr(c/(HVl^)f sinh^ x(l + z)]. (45) 

In the denominator we again have iirdpll + z), where dp is the proper distance (Weinberg 1972, 
MM95). The relation between x a-nd 1 -|- 2; is now given by 



1 + z = (coshrjo - l)/(cosh(?7o - x) - 1), 



(46) 



where cosh?7o = (2/i7) — 1. Hence 



sinh X = 7^(7——) (2(1 - ^) + ^(1 + ^) ± (2 - - n + n(l + z)) . (47) 




To decide the right sign we note that for z ^ one should have {cl{H\Jl — sinh x = cz/H. 
This is fulfilled for the minus sign. Hence here, instead of equation (37), one obtains 



i Fbii + zd') _ 2(1 + n(i + zd') - (2 - + + 

V Fd'{i + zb) ~ 2(1 + n(i + zb) - (2 - - n + n(i + zb) ' 



(1 + Zd,)/(1 + Zb) = r 



d'b ■ 



(48) 



The system is solvable analytically for a given tJ^^^Q. and Fb/F^' (e.g., substituting (1 + 2;;^/) 
from the second equation into the first one, a fourth order algebraic equation is obtained for the 
unknown (1 + -^fe)-) Nevertheless, it is much simpler to solve it numerically. In order to illustrate 
the trend, we solve it for T^'b = 2.25, j = 0.6, (Fb/Fd') = 25, for various i7. The relevant Zb and 
Zdi are summarized in Table 6. As seen from Table 6, decreasing values of i7 increase the relevant 
redshifts, as expected. More specifically, the numerical examples calculated here show that for the 
smallest values of i7 ~ 0.2 currently thought to be acceptable, an increase in redshift by up to a 
factor 2 may occur over the corresponding i7 = 1 values. 



We have derived analytical expressions for the integral distribution of bursting or steady 
sources, as well as the mean redshifts, time dilations and the dispersions of these quantities over 
finite fiux bandwidths, valid over the entire range of fiuxes in a spatially fiat Friedmann model. 
We have also evaluated the effects associated with i7 < 1 and with the inclusion of K-corrections. 
These expressions are particularly useful for comparing against the numerical results log N — log F 
fits, e.g. GRB catalogues such as the BATSE and PVO samples, and allow one to derive relatively 
simple and quick results without lengthy computations. 

We have used our analytic expressions for the average redshifts in specific peak fiux bands 
together with log N — log F fits to the BATSE and PVO data on GRB to derive new information, 
based on measurements of time dilation effects. While the detection of such cosmological 
time dilation signals is currently the subject of debate (Norris, et al. , 1995, Fenimore, et al. 
, 1995, Mitrofanov, et al. , 1994, Mitrofanov, et al. , 1995), there is no doubt that under the 
cosmological interpretation such effects should be incorporated into analyses of the integral 
distribution in order to obtain more reliable conclusions about the redshift, luminosity and density 
behavior. If a measured time dilation exists, several important problems can be addressed. First, 
the ambiguities of previous log N — log F fits (which allowed various choices of luminosities 



6. 



Discussion and Conclusions 
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and of the density evolution index _D, e.g. HMM95) are resolved. The inclusion of redshift 
information via a time dilation signal breaks the degeneracy of the models fits and allows a unique 
determination of £ and D. Second, the use of our analytical redshift expressions using parameters 
from log N — log F fits provide an estimate of the mean source redshifts in different peak fiux 
bands that are explicitly consistent with a particular logiV — log_F fit. Third, the same analytic 
expressions provide an explicit estimate of the dispersion (or essentially la error bars) associated 
with these mean redshifts. This dispersion is especially large for the examples of a power law 
luminosity function, which provide good log N — log F fits in a sense that are statistically 
indistinguishable from those of standard candles. 

We have also shown that a simple analytic approximation to the Fenimore, et al. , 1995 
energy stretching phenomenon of GRB time pulses implies that the measured time dilation factors 
between two arbitrary GRB fiux classes labeled by d" and h would scale with the redshift factor 
approximately as T^ih = {^td) I {^tb) = = [(1 + + zj^yf"^^- This is a particularly simple 

scaling that allows us to find theoretical estimates of the bright and dim mean redshifts of bright 
and dim GRBs, if a time dilation signal t^i}, is identifiable as being due to purely cosmological 
effect. Thus, for the value T^'b = 2.25 cited by Norris, et al. , 1995, we find that the redshift factor 

5 1 3 

ratio would the have to be r^'b = (1 + + ^b) = tJ^^ ~ 3.86, if energy stretching is included. 

The numerical examples above are for a stretching oc with a specific value (Fenimore, 
et al. , 1995) k = 0.4; we have in the text also provided general expressions for an arbitrary 
stretching index k. 

As specific examples, we have discussed the redshifts implied by reported values of the time 
dilation r^/fc = (Aid/)/(A4) ~ 2.25, 1.75, 1.35 (Norris, et al. , 1995, Norris, 1995b, Horack, et al. 
, 1995b, Fenimore, 1995b), under the assumption that these real and entirely cosmological, i.e. 
not contaminated by intrinsic effects (e.g. internal source physics which might mimic such a 
signature, etc.). Approximate standard candle (SC) results were obtained used the simple analytic 
equations 37, and more accurate results were obtained using the detailed analytic averaging over 
fiux bands corresponding to bright and dim+dimmest bursts. The latter are detailed in Table 1 
and Figure 4, from which one can find the relation between a particular redshift factor ratio r^t}, 
and the density evolution factor D (corresponding to a luminosity which gives consistency with 
the log A — log_F constraints). Table 1 and Figure 4 show that for t^i}, = 2.25, r^t}, ~ 2.25 (no 
stretching) occurs for D c:^ 4, where ^ ~ 0.5 and IJi ~ 2.3, while r^if, = 3.86 (stretching) occurs 
for _D ~ 5, where ^ ~ 1.6 and IJi ~ 10.0. For T^if, = 1.75 and no stretching one infers _D ~ 3 and 
'zb ~ 0.25, 'zj' ~ 1.25, while with stretching (r^'fe = 2.54) one gets _D ~ 4 and ^ ~ 0.5, 'zj' ~ 1.5. 
For Td'b = 1-35 and no stretching the implied value is _D ~ 1, ^ ~ 0.1 and IJ' — 0.7, while with 
stretching (r^'fe = 1.65) it is _D ~ 2.5, ^ ~ 0.25 and 1^7 ~ 1.15. 

The same three values of a putative time dilation Td^f, = (At(ii)/(Atb) ~ 2.25, 1.75, 1.35 
were also used for comparing with the band-averaged redshift factor ratios in the power-law (PL) 
luminosity function case, detailed in Table 3 and Figure 5. From these one sees that the average 
redshifts in the PL luminosity function case are not significantly different from those obtained in 
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the SC case case. For T^'b = 2.25 without stretching a value of r^'b — 2.25 would occur for _D ~ 4, 
where ^ ~ 0.9 and IJ' — 3.0; while with with stretching a value of r^'b = 3.86 occurs for _D ~ 5, 
where ^ ~ 2.0 and IJ' — 10.0. For T^'b = 1-75 (no stretching) we have D ~ 3.5, ^ ~ 0.6 and 
1^ ~ 1.8; with stretching (rd^b = 2.54) we have D ~ 4.5 and ^ ~ 1.4 and 1^ ~ 5.3. For T^^b = 1-35 
(no stretching) we have _D ~ 2, ^ ~ 0.3 and IJf ~ 0.85; with stretching (rd^b = 1-65) we have 
D ~ 3.4 and lb — 0.6 and IJf ~ 1.7. 

For the standard candle case without the energy stretching of Fenimore, et al. , 1995, our 
theoretical model fits to the data give mean redshifts of bright and dim sources. We find from 
our SC/no stretch fits that for T^'b = 2.25, 1.75, 1.35 the density evolution index required (see 
equation (1)) must be _D ~ 4, 3 and 1. With energy stretching, they would be _D ~ 5, 4 and 2.5. 
The large D for T^tb = 2.25 is in qualitative agreement with the results of Fenimore & Bloom, 
1995, Fenimore, et al. , 1995, Horack, et al. , 1995b in independent and complementary analyses. 
In the power law (PL) luminosity function, where we used a ratio of brightest to faintest intrinsic 
luminosity K = 10^, for the same three values of T^ib the no stretch fits are D ^ A, 3.5 and 2, 
while with stretching they are _D ~ 5, 4.5 and 3.4. The source luminosities corresponding to 
these density evolution exponents D (compatible with the log N — log F constraints) are given 
in HMM95. Here _D > 3 means more sources at large redshifts and _D < 3 means more at small 
redshifts, since in our proper density notation _D = 3 corresponds to comoving constant density. A 
confirmation of a clear cosmological time dilation would be needed to obtain stronger constraints 
on the density evolution. 

We note that the redshift values and evolution indices quoted are based on a statistical 
average over specific fiux ranges, giving mean redshift values which are systematically larger than 
the values (41) based on the simple analytic estimates of equations (37) for the centroid of the 
fiux range. Nevertheless, the dispersions of the average redshifts and their ratios are large, and 
therefore as a rough first approximation the values from (41) are compatible with those of Table 1 
for the SC case. For instance, subtracting a la dispersion from the mean value of r^tb for D = 5 
gives an r^^b comparable to the mean r^^b value for D = 4.5, corresponding to mean dim redshifts 
IJf ~ 5.5, if Td'b = 2.25. This is near the redshift of the most distant quasars known, and thus even 
for this most extreme time dilation in the SC case there need not be incompatibility between the 
Norris, et al. , 1995 signal and the earliest galaxy formation redshifts. 

The tendencies outlined above are stronger in the case of a power law (PL) luminosity 
distribution. While the mean values obtained are comparable, the redshift dispersion is 
understandably larger than in the SC case. For a relatively large ratio of maximum to minimum 
luminosities K = 10^, which still give good logiV — log_F fits, the PL dim and bright mean 
redshifts without or with energy stretching are somewhat larger than those in the corresponding 
SC cases given above. In the PL case, however, while again the mean dim redshifts are very large 
for a dilation signal of 2.25, whether stretching is included or not, the dispersion is so large that a 
dim redshift as low as z^i ~ 1 is within the la error bars of the mean values. 
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The smaller time dilation signals T^'b (e.g. 1.75 or 1.35) of course bring the mean redshift 
values closer to z^' ^ 1. There is also the additional possibility that the purely cosmological part 
of a time dilation signal be even smaller than what is directly measured, in that some specific 
models already imply a fairly generic source-intrinsic time dilation signal. As an example, the 
dissipative fireball external shock spectra calculated by Meszaros , Rees and Papathanassiou, 1994 
show that shorter duration bursts (higher T) also have higher intrinsic fiuences (because from 
general physical arguments one finds that the radiative efficiency generally increases with F), e.g. 
as illustrated in their Figure 4 of that paper. This mimics a cosmological time dilation, it is not 
strongly model dependent, and is independent of redshift (see also Brainerd, 1994, Yi, 1994). 

The specific numbers discussed above are for an i7 = 1, A = 0, K-correctionless model. For 
an i7 < 1 universe, or in cases where a K-correction is necessary, the mean redshifts are generally 
larger than in the i7 = 1 uncorrected case, but the error bars remain large enough that even time 
dilations ~ 2 can remain compatible in the PF case with modest dim source redshifts. 

In conclusion, we have derived analytic expressions for the integral number of sources per 
unit fiux in an i7 = 1, A = Friedmann model with density evolution and either standard candle 
or power law luminosity functions, of general applicability to steady or bursting sources. We 
have shown how the density evolution affects the luminosities inferred for the sources, and have 
discussed the effects of a power law luminosity function on log A — log_F distributions. Using 
these expressions we have also derived the mean redshift and time dilation factor distributions 
over given finite fiux bands, and the dispersions associated with them. We have shown that a 
time dilation signal of 2.25 (Norris, et al. , 1995), if purely cosmological, would imply values for 
the redshift of the dimmest bursts which are very large especially if one includes pulse energy 
stretching (Fenimore, et al. , 1995). However the redshift dispersion expected is so large that 
the redshifts remain statistically compatible with conventional ideas about the epoch of earliest 
galaxy formation. For smaller values of the time dilation signal, e.g. 1.75 or 1.35, even the mean 
values of the redshifts are in the conventional range. More generally, we have discussed the effects 
of density evolution and luminosity function on the redshift ranges inferred from an observed 
fiux distribution with various amounts of intrinsic pulse energy stretching for arbitrary values 
of an observed cosmological time dilation, which may be useful in the interpretation of future 
experiments. 

Acknowledgments: We are grateful to M.J. Rees and E. Fenimore for stimulating discussions, 
F Horvath and V. Karas and E. Fenimore for useful interactions, and J. Horack for detailed 
comments on the manuscript. This research has been supported by NASA NAG5-2857. P.M. is 
grateful for the hospitality of the Institute for Theoretical Physics, University of California, and 
support through NSF Grant PHY94-07194. 



7. Appendix: Details of equations (18-19) 
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One has 

n\ ( ''^ \ ra(ra — l)...(ra — A; + 1) 



kj \n — k J kl[n — ky. kl 
First, we prove the following relation (here ra > 0): 



, ra, A; = 0, 1, 2, k < n. 



E(-i)Mr)=o- (A.i) 



k=o 



Proof: If ra = 2a; — 1, where x is natural, one has ( — 1) = —( — 1)" for any k < n. Then 



ra \ / ra 
A; / ' \n — k 



occur with opposite signs, and they zero give in (A.l). We have x such pairs in the sum of (A.l); 
hence, here (A.l) holds. If ra = 2x, where x is natural, then one may write 



2x 

E(-i)* 



, L / 2x \ / 2x \ / 2x \ / 2x \ / 2x \ / 2x \ / 2x 

k=0 



k J \ J \ 2 J \2x J \ 1 J V 3 y V2a; - 1 



We have (x + 1) positive terms (the first and the last terms equal to 1), and x negative terms. We 
will use the following relation: 

ra\ /ra — 1\ /ra— 1\ 

k) " [ k ; ^ U- J ' 

where ra > and A; > 0. Hence 

|M)'(t)-.E[(1;Xr;)l-£[e:;)^(r; 

(The two factors of unity were written down separately; note the sum indices.) After some 
arrangement it follows that 

p-.'(t)-.E(ir)--t(:::;)-.E(-ir)-i:(-^>«. 

k=0 ^ ' m = l ^ ' m = l ^ ' m = l ^ ' r=0 ^ ' 

which completes the proof of (A.l). 

Using (A.l) it follows that (for ra > 2) 

r\i — J. r\i — J. ?Ti — U 



and (0 < s < ra - 2): 



^ (_l)'=fc(fc _ X)...{k -s)r)=n{n- l)...(ra - .) E U - I - 1 

k=s+l ^ ' fc=s+l ^ 
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n — s — 1 / -1 

n — s — \ 



m=u 



n(n-l)...(n-s) ^ (-l)"+^+M =0. (A.3) 



We define 8-2D = M; M = 0,1,2, 3, and we have (0 < A; < M): 

3(-l)'= fM' 



and thus (if M > 1): 

M 

Y,{k + ^)ak = 0. 

k=0 

In what follows we will take M > 2. The cases M = and M = 1 will be discussed later. One 
may write 

M k M -, 



fc=0 ^ fc=0 ^ " 

0' 



0,' 



2\ /2\ 1 /2\ 1 



+... 

J') '"TTi; + + '"(IT 



(Ik - 



yf 



a-M - \ ttM— h...+ (-l) ,^ «M- 



This "Christmas-tree" arrangement allows one to proceed as follows. First, consider the first terms 



of the horizontal lines. Their sum is 
M . , X M 



tTn \ ^ y tTn Jo Jo 



k=0 ^ ' k=0 

3 6 3 6 



Ad. 



M + 1 M + 2 M + 3 (M + l)(M + 2)(M + 3) 
Second, consider the second terms of the horizontal lines. Their sum is 

M ^ ^ N M M 

5*1 = - ^ f j ttk = - "^kuk = - ^(A; + 3 - 3)ak = SAd- 
k=i ^ ^ k=o k=o 
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(The lower index may be changed from 1 to zero, because we add only a vanishing term.) Similarly, 
we obtain for the third terms 

^2 = E (2) «^ = E '-^-^ = E = f:(-2)ka, = 6A0. 

k=2 ^ ^ k=2 ' k=0 k=0 

Consider now the s-th terms (s > 2). Their sum is: 

, ^ ( k \ i^ fc(fc-l)-(fc + 3-.)(fc + 2-.) 

, ^ - + 3 - s){h + 3 - s - 1) 
(-1) — «fc 



fc=0 



(.-1)! 



(-l)(-.-l) ^ , fc(fc-l)-(fc + 3-.) . + 1 

fc=o ^ ' 

^From this the following recurrent relation follows: 



m + 2 

Sm = Sm-x\ m = 1,2,...,(M- 1); 5o = A^. 

m 

This may be changed into the direct relation as follows: 

Sm = Sm-Xi Sm-\ = ':Sm-2-, ■■■ 5*2 = '^Si,Si = —Sq = SAjJ. 

m m — 1 2 1 

Hence 

S„ = "" + ^»"+" ,4,. (A.4) 

This relation holds for any < m < M . 

This "Christmas-tree" proof was done for M > 2. If M = 0, = 1 and 5o = 1. If M = 1, it 
is easy to show that = 1/4, Sm = A£)(m + l)(m + 2)/2, m = 0, 1. Hence (A.4) holds also for 
M = and M = 1. 

Using (A.4) one may write 

tUt t^k I- ^4^^ ^(k + l)ik + 2) , , , , , , 

E «fc n ^ f4u = / n ^ f4^3 E on ^ ,4U = ^D^Mit). (A.5) 



(1 + ^^)^ to (1 + ^')' Ml + t' f to 2(1 + 

To do the concrete integration one may proceed as follows (M = 0, 1, 2, ...). We write: 

t^dt 



Mi) 



(1 + ^4) 



3 ■ 



f5 33 

Mi) = j^^^s + T^v^o(0. 
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The function ipo(t) is a tabulated integral. Consider now M > 2, and assume that we know 
tpM-iit) and tJjM-2i't)- Then one may write: 



Band, D., 1994, ApJ, 432, L23 

Brainerd, J.J., 1994, ApJ, 428, LI 

Cohen, E. and Piran, T., ApJ, in press 

Dermer CD., 1992, Phys. Rev. Letters, 68, 17 

Emslie A.G., Horack J.M., 1994, ApJ, 435, 16 

Fenimore, E.E., et al. , 1993, Nature, 366, 40 

Fenimore, E.E., et al. , 1995, preprint 

Fenimore, E.E., 1995, private communication 

Fenimore, E.E. & Bloom, J.S., 1995, ApJ, 448, LlOl 

Fishman G.J. et al. , 1994, ApJS, 92, 229 (BATSE IB catalogue) 

Hakkila, J., et al. , 1994, ApJ, 422, 659. 

Hakkila, J., et al. , 1995, Proc. ESLAB Symp., in press 

Horack J.M., Emslie A.G. & Meegan C.A., 1994, ApJ, 426, L5 

Horack J.M., Emslie A.G.& Hartmann D., 1995, ApJ, 447, 474 

Horack J.M., Mallozzi R.S. k Koshut T.M., 1995, preprint 

Horvath, L, Meszaros , P. and Meszaros , A., 1995, ApJ, in press (astro-ph/9507053) 

Loredo, T. & Wasserman, L, 1992, in Compton Gamma-ray Observatory, eds. Friedlander, M., 
et.al, (AIP Conf.Proc. 280, Washington, D.C.), 749 
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Table 1 



D 




(1 + ^6) 




(1 + Zd') 




rd'b 




1.0 


0.2 


1.15 


0.036 


1.59 


0.09 


1.39 


0.09 


1.5 


0.3 


1.18 


0.046 


1.74 


0.11 


1.47 


0.11 


2.0 


0.45 


1.23 


0.057 


1.94 


0.14 


1.58 


0.14 


2.5 


0.6 


1.26 


0.066 


2.11 


0.17 


1.67 


0.16 


3.0 


0.8 


1.31 


0.078 


2.32 


0.21 


1.78 


0.19 


3.5 


1.1 


1.36 


0.093 


2.61 


0.26 


1.91 


0.23 


4.0 


2.0 


1.50 


0.132 


3.36 


0.40 


2.23 


0.33 


4.5 


7.0 


2.02 


0.287 


6.55 


1.04 


3.24 


0.69 


5.0 


15.0 


2.63 


0.474 


10.90 


1.96 


4.14 


1.05 



Table 1 Caption: Standard Candle results based on equations (38, 39) of the mean {^-\-z) values 
and their ratio and dispersions for different D, where the different Co (in units photons/ (err? s)) 
used for different D are taken from the fits of HMM95. 
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Table 2 



D 




(1 + ^6) 








Td'b 




4.5 


5.0 


1.84 


0.23 


5.38 


0.80 


2.92 


0.57 


4.5 


6.0 


1.93 


0.26 


5.98 


0.92 


3.09 


0.63 


4.5 


7.0 


2.02 


0.29 


6.55 


1.04 


3.24 


0.69 


4.5 


8.0 


2.10 


0.31 


7.11 


1.16 


3.38 


0.75 


4.5 


9.0 


2.18 


0.34 


7.65 


1.28 


3.51 


0.80 



Table 2 Caption: Standard Candle results showing the effects on the mean redshifts of allowing 
the value of Co (ph cm~^ ^~^) to vary within ^ Isa of its best fit value for fixed D. 
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Table 3 



D 




1 (1 + ^6) 


A(l + ^6) 




1 A(l + ^d,) 


rd'b 


Ar^/fc 


2.00 


0.05 


1.33 


0.22 


1.82 


0.63 


1.37 


0.52 


2.50 


0.08 


1.42 


0.28 


2.07 


0.86 


1.46 


0.67 


3.00 


0.10 


1.49 


0.33 


2.33 


1.12 


1.57 


0.83 


3.50 


0.15 


1.62 


0.42 


2.83 


1.63 


1.75 


1.11 


4.00 


0.30 


1.92 


0.66 


4.06 


2.99 


2.12 


1.72 


4.50 


0.50 


2.32 


0.98 


6.29 


5.34 


2.71 


2.57 


5.00 


0.80 


2.96 


1.46 


11.00 


9.63 


3.70 


3.73 



Table 3 Caption: Power Law luminosity distribution results showing average redshifts 
calculated from equations (38, 39), notation as in Table 1. A ratio of K = Lm I J^m = 100 was 
taken everywhere, and values of Cm (ph cm~^ ^~^) for each D are taken from the best fits of 
HMM95. 
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Table 4 



D 


£„/(105^/.-2) 


(1 + ^6) 


A(l + ^6) 




A(l + ^dO 


rd'b 


Ar^/fc 


4.50 


0.20 


1.84 


0.56 


4.22 


2.82 


2.29 


1.69 


4.50 


0.40 


2.18 


0.85 


5.68 


4.56 


2.60 


2.32 


4.50 


0.60 


2.44 


1.09 


6.86 


6.08 


2.81 


2.79 


4.50 


0.80 


2.66 


1.31 


7.91 


7.48 


2.97 


3.17 


4.50 


1.00 


2.86 


1.51 


8.87 


8.81 


3.11 


3.49 



Table 4 Caption: Mean redshifts for the power law luminosity function; D and K are fixed; 
K = Lm I J^m = 100 and Cm (ph cm~^ ^~^) varied. 
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Table 5 



D 




K 


(1 + ^6) 




(1 + Zd') 




rd'b 




4.5 


0.5 


4 


1.37 


0.12 


2.61 


0.47 


1.90 


0.38 


4.5 


0.5 


64 


2.10 


0.74 


5.51 


3.95 


2.63 


2.10 


4.5 


0.5 


324 


3.12 


2.04 


8.88 


11.4 


2.85 


4.10 


4.5 


0.5 


1024 


4.29 


4.14 


12.2 


23.0 


2.84 


6.03 


4.5 


0.5 


2500 


5.51 


7.11 


15.3 


39.1 


2.77 


7.94 



Table 5 Caption: Mean redshift in the case of a power law distribution. Cm (ph cm ^ s 
and D are fixed, and K = Cm I J^m varies. 
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Table 6 

n 1.0 0.8 0.6 0.4 0.2 

Zb 0.6 0.7 0.8 1.1 2.0 

5.2 5.5 6.0 7.0 10.6 



Table 6 Caption: Values of the mean resdhifts derived for various values of^l. 



- 32 - 



Figure Captions 

Figure 1. Standard candle (SC) theoretical integral peak flux distributions for values 
D = 2, 3, 4, 5, 6, 7 in increasing order upwards. Two sets of curves are shown, for two values of 
the luminosity Cg- The bottom set has a luminosity Cg = 4.5 X 10^^h~^ s~^, the upper one has 
= 3 X s-^ The lower set has been artiflcally moved down by one unit in log N so as 

not to overlap. For comparison, the 2B data set is also shown twice, again the lower curve down by 
one unit, to illustrate qualitatively that different choices of £o provide approximate fls to the same 
data set corresponding to different evolution indices D. The higher luminosities flt with higher D. 
In this case the lower curves flt approximately to _D = 2, the upper ones with D = 4 (and similar 
flts can be found for _D = 3. More accurate flts are discussed in HMM95, see also Table I. 

Figure 2. Power law luminosity function (PL) theoretical integral peak flux distribution 
for D = 2,3,4,5,6,7 increasing upward, for [i = 15/8. Two sets of curves are shown, both 
for K = Cm / = 10^. The lower set (which has been offset downwards by one unit) is for 
(CM,^m) = (5 X 10^'' s~^,5 X 10^^/i~^ s~^) and suggests a flt for D = 2, while the upper set 
has (CM,^m) = (2.5 X lO^^/i"^ s~^,2.5 X lO^'^/i"^ s~^) and suggests a flt to _D = 4, when compared 
to the 2B data set. Again, higher luminosities provide qualitative flts to higher evolution indices 
D. For more detailed flts, see HMM95 and Table 3). 

Figure 3. A comparison of theoretical integral peak flux distributions for standard candle 
(SC, upper curves) and power law (PL, lower curves) luminosity functions. The curves in each 
set are, in upward order, for density evolution indices D = 2,3,4,5,6,7. Arbitrary values of 
Co = W^^h-^ s-i for SC and {CM,Cm,K) = (W^^h-^ s'^W^^h-^ s'^W^) for PL cases were 
chosen. Both SC and PL show (for D < 4) two regimes, N(> F) oc F"-^/^ at high F and F 
constant at low F. The large ratio K = Cm I J^m in the PL case was chosen to illustrate the third 
asymptotic regime in this case at intermediate _F, namely N{> F) oc _F^~'^, where fi is the index 
of the power law luminosity function. 

Figure 4- Standard candle luminosity function: Mean (1 + z) for dim+dimmest (bottom 
curve) and bright (upper curve) bursts, and their ratio r^'b (middle curve), as a function of the 
density evolution parameter D. Values are computed for integer and half-integer values of D, but 
are drawn slightly offset in the flgure so the error bars can be distinguished. For a time dilation 

5 /3 

ratio Tf^/ft one has a redshift factor ratio r^/f, = t^i}, {r^t}, = tJ^^ ) if one neglects (includes) energy 
stretching of time proflles. Results based on equations (38, 39), as in Table 1, where the different 
Co used for different D are taken from the SC flts of HMM95. 

Figure 5. Power Law luminosity distribution: Mean (1 + 2;) for dim+dimmest (bottom curve) 
and bright (upper curve) bursts, and their ratio r^'b (middle curve), as a function of the density 
evolution parameter D. Values are computed for integer and half-integer values of D, but are 
drawn slightly offset in the flgure so the error bars can be distinguished. For a time dilation 

5 /3 

ratio Td'b one has a redshift factor ratio r^^b = Td'b (I'd'b = Tjb ) neglects (includes) energy 

stretching of time proflles. The results are calculated from equations (38, 39), as in Table 3, using 
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a ratio of maximum/minimum luminosities K = 100 and values of for each D taken from the 
PL fits of HMM95. 
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